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Abstract 

Elementary particle scatterings and decays in presence of a background magnetic field are very 
common in physics, specially after the observation that the core of the neutron stars can sustain a 
magnetic field of the order of 10^^ G. The important point about these calculations is that they are 
done in a background of a gauge field and as a result the calculations are prone to gauge arbitrariness. 
In this work we will investigate how this gauge arbitrariness is eradicated in processes where the initial 
and final particles taking part in the interactions are electrically neutral. Some comments on those 
processes where the initial or final state consists of electrically charged particles is presented at the 
end of the article. 

1 Introduction 

Calculations of elementary particle decays and scattering cross-sections in presence of a background 
magnetic field are commonly found in literature m21iniEl- These calculations became more important 
after it was understood that the neutron star cores can sustain magnetic fields of the order of 10^'^ G or 
more. In presence of such strong magnetic fields all the particles which have a magnetic moment are 
bound to get affected and consequently their properties like self-energy, decay rates or scattering cross- 
sections are modified. The interesting feature of a background magnetic field is that it not only modifies 
properties of particles with magnetic moments but can also affect the properties of particles which do 
not have any magnetic moment. The obvious question is, how is it possible? To give a satisfactory 
answer to this question we take the example of standard model neutrinos and their self-energy. In the 
standard model, neutrinos have no electric charge hence they do not have any direct coupling to photons 
in any renormalizable quantum field theory. The standard Dirac contribution to the magnetic moment, 
which comes from the vector coupling of a fermion to the photon, is therefore absent for the neutrino. 
In the standard model of electroweak interactions, the neutrinos cannot have any anomalous magnetic 
moment either. The reason is: anomalous magnetic moment comes from chirality-flipping interactions 
'ipo'tii^i-'F^'^ , and neutrinos cannot have such interactions because there are no right-chiral neutrinos in 
the standard model. Consequently the standard model neutrinos cannot interact with the background 
magnetic field. If one looks at the Feynman diagrams contributing to the neutrino self-energy in Fig. ^ 
one will see that the self-energy diagrams contain charged lepton, £, propagators and W~ propagator. 
Although the neutrinos in this case have no coupling with the background fields the virtual particles in 
the loops can have couplings with the background field and so the result of the self-energy calculations 
arc bound to be affected by the background fields. 
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Figure 1: One-loop diagrams for neutrino self-energy. 




Figure 2: One-loop diagram for vacuum polarization of the photon. 

As a second example we take the case of the vacuum polarization of the photon. To one-loop the 
Feynman diagram is shown in Fig. [21 In this case also the virtual particles propagating in the loops must 
be charged particles and as a result the magnetic field will affect the expression of the vacuum polarization 
result. There can be various other cases, as an example the four photon interactions in QED, in which 
the magnetic field can affect properties of electrically neutral particles. 

The calculations of the quantities like the neutrino self-energy and the vacuum polarization of the 
photons in presence of a background magnetic field is similar to the calculations of them in absence of 
the same, except that now one has to use the modified two-point functions of the charged particles. We 
will discuss about the charged fermion two-point functions in presence of a magnetic field in sectional 

Before going into further details about how to calculate various elementary particle processes in 
presence of a background uniform magnetic field, we must be careful about the fact that we are calculating 
physical quantities in presence of a gauge field. The background field must be specified by some suitable 
choice of gauge. Then the question is whether the quantities which we calculate are in the end gauge 
invariant? If the answer is no than there is no utility of doing all those calculations because they will be 
arbitrary to the gauge which one uses. In this article we will try to understand how the gauge arbitrariness 
of the various calculations mentioned above can be tackled. We will mainly focus on all those elementary 
particle processes whose Feynman diagrams do not have any charged particle external legs. Processes 
like Compton scattering, pair annihilation and so on will be dealt briefly in the penultimate section. 

The article is organized in the following way. In section |21 we will set the preliminaries required to 
follow the following sections and doing so we will also set the mathematical notations which we will 
follow. In section |21 we will discuss about the form of the Schwinger phase accompanying the charged 
fermion two-point function in presence of a magnetic field and try to understand its properties. The 
next section will deal with the effects of the phase factor on loop calculations. It will be discussed in 
which circumstances we can neglect the phase factor for an actual calculation and when this cannot be 
done. In sectionElwe show the Dirac equation solutions in a magnetic field background and a pure gauge 
background and comment on their gauge transformation properties and the difficulties we come across 
to prove background gauge invariance of those processes which involve charged particles in the external 
lines of the Feynman diagrams. 
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2 Conventions and notations 



In presence of a background magnetic field we can decompose the photon field as follows: 

(1) 

where A'^{x) is the dynamical photon field and A'^{x) is the classical background field which gives rise 
to the magnetic field. If the uniform background classical magnetic field is called B then we must have: 

B^VxABix), (2) 

where A'^{x) — {0, Ab{x)). In presence of the background magnetic field we can also write the field 
strength tensor as: 

F'^^x)^F^''{x)+F^\ (3) 

where F^^^ix) = - d''A'^{x) and F^^ = d'A^{x) - d^A^{x) is a constant as given in Eq. 

The QED Lagrangian can be written as: 

£ = V^(i7M^^ - - i^^'^M- > (4) 

where D^^ = 9^ + ieqA'^ is the covariant derivative of the fermion fields. Here e is the magnitude of the 
electronic charge and q designates the sign, for electrons q = —1 and for positrons q = +1. In this article 
we will only be talking about those cases where |g| = 1. The gauge invariant Lagrangian can also be 
written as: 

£ - ^ [7^n^ - m] - eq-^j^M^ " Ipf^" F^. , (5) 

where 11^ — id^ — egAg is the kinetic momentum of the fcrmions in presence of the background field. 
The commutation relations of interest are: 

[x^'.n"] = ^ig^'" , and [W, H""] = -leqF^'' , (6) 

where, 

3^'' = diag(l, -1,-1,-1). (7) 

Now we can talk about the gauge invariance of the Lagrangian given in Eq. Q . We have to be particularly 
careful since there are two kind of gauge invariances here. The first kind is about the gauge invariance 
of the Lagrangian under a gauge transformation of the dynamical photon field A'^ —^ A'^ + d^uj where 
u! is some well behaved function of space and time. Under this kind of a gauge transformation the 
fermion fields will transform like ■0 ■06"''"''^. But this does not exhaust all the gauge transformation 
possibilities of the Lagrangian given in Eq. ((nj, we can still have a gauge transformation of the background 
gauge field Aq as Aq A'^ + d^X where A is a well behaved function of space-time coordinates, which 
will leave the Lagrangian in Eq. ((SJ invariant provided the fields transforms like tpe~'^^'^^. 
The Euler-Lagrange equation of ip is: 

{fl - m)^/j{x) = , (8) 

where fl = jfj.H'^- Generally the inverse of the kinetic part of the Lagrangian (for Eq. the quantity 
sandwiched between the fields in the first term in the right hand side of the equation) gives the two-point 
Greens function of the theory, which are the building blocks of subsequent developments. In the present 
circumstance, the two-point functions are also defined as: 

{fl-m)SB{x,x')^d^{x--x'), (9) 
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which can also be written in a matrix form: 



Sb = ^ , (10) 

/i — m 

where Sb is a matrix, in fact the inverse matrix offl—m, and the two-point Greens function is (a;'|S'B|a;) = 
Sb{x,x'). So from the definition of the two-point function given in Eq. we see that the fermion two- 
point functions are defined not in a gauge invariant way but has exphcit dependence on the background 
gauge field. We will talk about this gauge dependence of the two-point function in the next section. 
In the Feynman diagrams appearing in Fig.^and Fig.|21 the photons represent the dynamical photons 
appearing in the Lagrangian given in Eq. . The background field A'^ is tacitly taken into care by the 
fermion propagators or other charged particle propagators appearing inside the loops. Another interesting 
point which must be noted is that although the Lagrangian as written down in Eq. |SJ| is Lorentz invariant 
but it contains A'^ and A'^ satisfies Eq. Q with the condition that A'^{x) = 0. The condition as given 
in Eq. |(5J) is not Lorentz invariant i.e., a magnetic field is not a Lorentz invariant quantity. So if we 
demand that we will be working in presence of a uniform classical background magnetic field then we are 
restricting the Lorentz covariance of the theory. 

Lastly we fix the notations and conventions which we will be following in this article. In future the 
magnitude of the uniform background magnetic field will always be denoted by the symbol B and we will 
take it to be pointed along the z-direction of the coordinate axis. The magnetic field vector pointing in 
the z-direction will be given by 

B^zB, (11) 

where z is the unit vector along the z-axis. 

As Lorentz covariance is restricted and the magnetic field chooses a particular direction in space, the 
4-vector structure breaks down into a perpendicular part and a parallel part. If is a 4-vector, then 

a|[ = (a°,0,0,a3), (12) 
= (0,a\a2,0), (13) 



such that 



Also in our convention 



where 



Terms as an and a , stands for: 



such that 



(14) 



df^i^^glu + a^u, (15) 



gl = (1,0,0,-1), (16) 

gi^. = (0,-1,-1,0). (17) 



a| = 4,a''a% (18) 

- (aY-(«^)^ (19) 

«i - -.9^.fl^fl% (20) 

= iar + {ar, (21) 



= ajj - ai . (22) 
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Also in our case, 



F^^ = -Fl"- = B . (23) 

With this amount of prchniinary understanding about the fermion two-point functions in a background 
magnetic field and the notations which we will follow we move to the next section which deals with the 
explicit form of the two-point fermionic Greens functions. 



3 The phase factor appearing in the two-point Greens function 
of a charged fermion 

In a much celebrated paper Schwinger derived what will be the form of the fermionic two-point function in 
presence of constant electromagnetic fields [HIIHIIZI- In the present article we will only study the behaviour 
of the two-point function in presence of a pure magnetic field. Later on it was possible to find a curious 
momentum space description of the Schwinger's two-point function in presence of a magnetic field. It can 
be noted that the above two-point function is not a free field propagator, it contains interactions with the 
background magnetic field. This fact makes the two-point function to be translationally non-invariant. 
The reason why it becomes translationally non-invariant is closely linked to its gauge transformation 
properties which we will discuss shortly afterwards in this section. As translational invariance is lost we 
cannot ideally Fourier transform S-q(x^x') into something analogous to S-q{p)^ but the actual two-point 
function contains a multiplicative factor which is translation invariant and has all the good properties 
that a two-point function should have. 

Schwinger's calculation of the fermion two-point function relies on solving the operator formula given 
in Eq. H10|l and then finding out the matrix element (x'jS'Bla:^) = Sb{x,x'). The two-point function can 
be expressed as |71 151 li^ 1111) : 



'4„ 



iSb{x,x') = k{x,x') J -0^e-'^P<^--'hSBip) , (24) 

here x stands for the coordinate 4- vector as usual. S'b(p) is expressed as an integral over a variable s, 
usually (though confusingly) called the 'proper time': 

/•oo 

iSb{p) = / dse*(P'^)G(p,s). (25) 







The quantities ^{p, s) and G(p, s) can be written in the following way : 



$(p,s) = is I pit - pl-m^] ~es, (26) 

' " eqBs ' 



Gip,s) 



cos{eqBs) \ " cos(ei3s) 
= (1 + iS^ tan(egBs))(^|| + m) + sec^{eqBs)^± . (27) 



In the above expressions = ^7 7 , where the gamma matrices are taken to be in the standard Dirac- 
PauH representation. — g^uj'^p'^ and — gj^^^^p" while the symbols p| and p\ are explained in the 
last section, e is an infinitesimal positive quantity introduced for the convergence of the integrals. For 
convenience henceforth we will call the expression of the two-point function as given in Eq. H24|l as the 
Schwinger two-point function. 

In a typical loop diagram, one therefore will have to perform not only integrations over the loop 
momenta, but also over the proper time variables. From the above equations we see that iSB(j>) is 
manifestly translation invariant, and it has another interesting property. As B — > 0, ^{p, s) — > is(p| — 
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— m^) — es = isijP' — w?) — es and G{p^ s) — > + m) H-^_l = |5 + m. Therefore when _B ^ we have: 

iSB{p) I dse'"(p'-™')-"(;5 + TO), (28) 







which is the normal fermion propagator in absence of any background magnetic field. More over it is 
seen that iSB(j>) consists of Bs and not the gauge fields, so «S'b(p) is not only translation invariant but 
gauge invariant also. This leads us to the conclusion that the gauge dependence of the two-point function 
iSB{x,x') must present in the function k{x,x'). 

Not going into any detailed description of k{x,x') we can simply understand its necessity in the two- 
point function. In presence of a background gauge field the gauge transformation property of the fields 
of the charged fermions are different at two different space-time points. Unless there is some factor in 
the two-point function which can connect these two fields at different space-time points with different 
gauge transformation properties, the calculations involving charged fermion two-point functions will not 
be manifestly gauge covariant. The gauge transformed fields comes with phase factors where the phase 
depends upon the space-time point where the gauge transformation is made. The fermionic fields at two 
different space-time points will therefore have two different phase factors. To make a connection between 
them k{x, x') must also be some form of a phase. Conventionally it is named the phase-factor. The phase 
factor, as calculated by Schwinger [S], is given as: 

k{x,x') — e-Kp {ieql{x, x')} (30) 

where 

i{x,x')= r 



K{0 + \F^\i- 
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(31) 



From Eq. (|31|) we notice that the phase factor breaks the translation invariance of the two-point function. 
For a constant background field we can always write the gauge field as 



^^(0 = -oi^re. + 9''A(e), (32) 



where A(^) is an arbitrary well behaved function and depends upon our choice of gauge. Using the above 
relation in conjunction with Eq. H31|l we can simplify the integration appearing in the phase factor as 



I{x, x) = / 



-iFrx; + a^A(e)] . (33) 

Using the constancy of the field strength tensor the above expression can be written as 

I{x, x') = ^x'^F^-x, + \{x) - \(x') . (34) 

From Eq. (|34|l we can immediately see if we set a; = x', in other words if we integrate over any closed 
contour in space-time I{x, x') vanishes. Thus I{x, x') connecting two points in space-time is independent 
of the path joining them, and as a result the phase factor of the Schwinger two-point function joining the 
points x' and x in Eq. (|3U|) is also path independent. 

Utilizing the path independence of the phase factor of the two-point function the general convention 
is to choose a straight line path connecting the two points x' and x. Points on this path are represented 

by 

= (1 - 0^'^ + C:^" , (35) 



6 



where the parameter ( ranges from to 1. Using Eq. and the straight hne path given in Eq. (|35|l 
one gets 



I{x,x') 



J x' 



(36) 



Using Eq. for the gauge field we can retrieve Eq. 

The form of Eq. (|3()|l when I{x, x') is as given in Eq. H36|l is very similar to a Wilson line but it is not 
a general Wilson line. A general Wilson line in presence of an U(l) field is defined as: 



U{x,x') 



exp 



leq I d^M^B 
p 



(37) 



where P specifies a particular path joining the points x and x' in space-time. From Eq. (|36l) we see that 
in this case I{x, x') is obtained by integrating Ag(^) only on a straight line path joining x' and x as given 
in Eq. (|35|l . For a general Wilson line as defined in Eq. (|37|l the gauge invariant Wilson loop will be: 



U(x, x) = exp 



(38) 



where 5 is a surface that spans the closed integration loop and dcr^i^ is an area element on this surface. 
As such U{x,x') will not be unity in general but as the integration defining I{x,x') in Eq. H36|l is only 
on a straight line path so in this case k{x, x) will always be unity. 

From Eq. (|34|l it is clear that the phase factor is dependent on the form of the function A(^), that 
is to say the fermion two-point function is dependent on the gauge in which the constant background 
magnetic field is specified. Suppose we are working in such a gauge that A(^) = 0, and then we make a 
gauge transformation of the background field as 



A^->A^ + a^A(e) 

then from Eq. (|24() it follows that the fermion two-point function will transform as 

iSB^x, x') ~* exp(jeA(x))iS'B(a;, x') exp(— ieA(a;')) . 



(39) 



(40) 



under the gauge transformation. As because we are working in presence of a background gauge field, the 
fields of the charged particles and their two-point functions both become background gauge dependent. 
This is the reason why the phase factor arises in the expression of the two-point function. 

In the B ^ {) limit iS-B{p) goes to the normal fermionic propagators in absence of any magnetic field, 
but what is the fate of k,{x, x') when B ^ 07 The answer is closely related to the way one choose Ag. 
There are many equivalent ways of writing as: 



AO=A^=A^=0, A^^-yB + b. 



(41) 



or 

A%=A^=A^^=0, Al=xB + c. (42) 

or 

A^=A| = 0, Al = ]^xB + d, = + (43) 

where 6, c, d and g are constants. In the above equations x, y are just coordinates and not 4- vectors. All of 
these above choices gives a magnetic field along the z-axis of the proper magnitude and the point to note 
about the above choices of the gauge fields is that all of them goes to zero as i? if 6 = c = d = g = 0. 
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When b = c — d — g — and i? ^ it is noted that iSB{p) approaches the free propagator form 
as it should be in absence of the gauge fields and k{x,x') also approaches unity, as is expected. But 
b = c = d = g = Ois not the most general gauge choice which produces a magnetic field along the 
^-direction. With the general choice of the gauges as is specified in Eq. H41|l . Eq. (|42|) or Eq. H43|l it is 
clear that k{x, x') does not approach unity as B ^ 0. If we take the choice of the background gauge field 
as given in Eq. H41I) then k{x, x') g^eqb{x-x ) ^ gQ -^^i^^Qn _B ^ in this case we will have, 

iSB{x,x')->e'''^''^^-^'^ I -^e-'P-("-"') ^(^ + "') , (44) 
J (27r)* — rn^ + le 

which is translation invariant but not the form which we expect when there is no background magnetic 
field present. 

A some what similar situation arises when we deal with pure gauge configurations. The Schwinger 
two-point function as given in Eq. H24|) is specific for the case of those gauges which gives a background 
magnetic field. In the general derivation by Schwinger W the form of the fermionic two-point function 
was derived for general constant F^'^ . If we choose a pure gauge field which does not give rise to any 
electric or magnetic field as 

= (45) 

where fc is a constant, then also we can find the Schwinger two-point function and it differs from the 
vacuum propagator only by the phase factor. Using Eq. (|35|l in this case we will have. 



iSB{x,x') = e(-<J'=(-°— / !^e-^M^-^')^yiLLp_ . (46) 



Jieqk{x°'x^-x"'x'J) f ^-'^P ^-ip-(x-x' ) ^(^ + ^71) 

27r^ — m? + ie 

This form of the propagator is not translation invariant and does not match our expected form. The 
exponential factor multiplieng the propagator of the free fermions is annoying. In the next section we 
will explicitly show how to tackle these problems. 

A brief summary of the properties about the Schwinger two-point function is presented below. The 
Schwinger two-point function, in general, is a function of A'^ and is defined for all gauge configurations 
including pure gauges. In the present section only its form in presence of a uniform magnetic field was 
presented. As a consequence of this fact when one takes the _B — s- limit of it, the free-fermion propagator 
is reproduced up to multiplicative phase, which designates a trivial pure background gauge configuration. 
In other words the B ~^ limit of the two-point function means a transition (not a gauge transformation) 
of iSB{x,x') where initially the gauge fields followed Eq. ((SJ with A'^{x) = and finally A'^{x) — dp{x) 
for some well behaved function p{x). The gauges which produces a magnetic field along the z-axis may 
contain constant terms as &, c, d, g and when _B = all the gauges in Eq. H41|) . Eq. H42|) and Eq. (|43|l 
becomes trivial pure gauge configurations. So in general the Schwinger two-point function as given in 
Eq. H24|) can be continuously transformed from the case where it is a function of the gauge fields which 
produces a magnetic field to the case where it is a function of pure gauge fields. This fact has an interesting 
outcome. It is known that in presence of an external magnetic field the transverse momenta p±^ do not 
represent gauge invariant degrees of freedom. But still the expression of the Schwinger two-point function 
contains p\ explicitly. The presence of p\ is there only because the two-point function is a function of 

(a;) and the domain of {x) which can give rise to a magnetic field can be continuously transformed 
into trivial pure gauge configurations. 



4 Effects of K,{x,x') on calculations 

As we have discussed previously in section |21 the phase factor in the Schwinger two-point function 
appears because the fermion two-point function attaches two points with different gauge transformation 
properties. Till now we have only talked about fermion two-point functions but one can also find out the 
Greens functions of charged scalars and vector particles in presence of a uniform background magnetic 
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field. It is seen that all of these two-point functions can be written as a product of a translationally 
invariant part and a background gauge dependent part. In fact all of the gauge dependent parts of the 
various two-point functions are same and are equal to k{x,x') 11 . The charged gauge boson two-point 
functions also depend upon other gauge parameters, like Feynman gauge parameter or Landau gauge 
parameter, but those are related to the dynamical gauge invariance. As we are not interested in an 
actual loop calculation we will not explicitly write down the forms of all the two-point functions here but 
for further discussions we will only utilize the fact that the translationally non-invariant part of all the 
two-point functions of charged particles are functionally equivalent to k{x,x'). So all the points wc will 
use to prove that a two-point fermion loop is independent of the choice of the background gauge will also 
apply for a loop like the one which is to the left in Fig. ^ 

In this section we will study various cases and try to understand what will be the effect of k{x,x') 
on loop calculations. To understand its importance we take an example of the one loop photon vacuum 
polarization in QED. Let P and Q be the two space-time coordinates where the photon line interacts with 
the virtual charged fermions. If we are interested in finding out the overall phase factor accompanying 
the vacuum polarization tensor then we will have to use the Schwinger two-point function for the charged 
fermions. The contribution from the phase factors k(Q, P) and k(P, Q) to the loop integral in coordinate 
space, which we denote as ^{P, Q) will be according to Eq. H30|l and Eq. (|34|l 

$(P,Q) - KiQ,P)KiP,Q) 

= expl^teq^[P,,F^''Q, + Q^F^''P,]^ , (47) 

which reduces to unity because of the antisymmetry of F^'^ . From Eq. (|47|) it is seen that the phase 
factor's contribution in the one loop calculation is trivial and obviously gauge invariant. The same 
analysis also holds for the diagram in the left of Fig. ^ 

Next we take a loop with three charged particle two-point functions connecting three space-time 
points P, Q, R. The overall phase contribution to the loop integral in the coordinate space can then be 
calculated using Eq. (|Sn|) and Eq. and is given by 

$(P,Q,i?) = k{Q,P)k{R,Q)k{P,R) 

= exp jiegi [P^F^^'Q, + Q^^F^'R, + R^F^'Pu] | • (48) 

As all the phase factors are of the same form the first point to notice is that the contribution from the 
function A(^) cancel out in the overall factor, showing that the contribution is explicitly gauge invariant. 

The next point which requires to be discussed is about the path independence of the phase factor. 
From section |31 we know that k(Q,P) is independent of the path which joins them, but here the path 
independence of k{Q, P) does not imply 

$(P, Q, R) = K(g, P)k(P, Q)h{P, P) = 1 , (49) 

or from Eq. 

/(P, P) = /(g, P) + /(P, Q) + I{P, R) = 0. (50) 

Instead of the above expectation we get a finite contribution from Eq. H48|) . As I{Q,P) consists of the 
product of the two end points instead of their difference, the different phase factors from the different 
paths connecting two intermediate points of the loop when multiplied does not reduce to unity. In general 
if, 

I{Q,P) = f{Q)-.f{P), (51) 

where /(P) is some well-behaved function of space and time, then Eq. H5U|) will hold and consequently 
Eq. H49() will be true. In this regard we can say that the factor which multiplied the the normal fermion 
propagator in absence of any magnetic field in Eq. 1I44|I is of the form as given in Eq. IjSlf) and consequently 
they will not pose any problems for actual loop calculations. 
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4.1 Similarity of $ with the Wilson loop 

In the previous discussions it was shown that the phase factor contribution to the one loop calculations of 
various cases are explicitly gauge invariant. In the case of the photon vacuum polarization it was shown 
that it contributes nothing for the phase factor. The contribution from the phase factors for cases where 
one has three or more than three vertices can be understood in another way. If we have a loop with three 
or more vertices then the overall phase will be c?^^^g(^) where L designates the path which joins all 
the vertices in straight lines. Then from generalized Stokes theorem we can write: 



where S is the area of the loop enclosed by the straight lines and dcTp^i, is the infinitesimal surface area 
in the — plane and F^'^ is the field strength tensor. In this case we notice that k{x,x') can be 
interpreted as the Wilson loop which is a gauge invariant quantity. 

For pure gauge configurations we have seen in Eq. H4t)|) that the fermion two-point function in presence 
of a pure gauge field comes with an unwanted phase. In the level of two-point functions this problem 
remains but if we look at loop integrals then this problem of the unwanted phase disappears. For the 
specific gauge choice in Eq. 14511 we see that /(P, Q) is of the form as given in Eq. (|51|l and consequently 
the overall phase factor will be trivial. But the result is not dependent on the particular form of the 
pure gauge chosen. From Eq. H52|) we can see in all those cases where one calculates some process which 
contains three or more than three vertices in presence of a background pure gauge field $ = 1 because 
for a pure gauge configuration F^'^ {S^) — 0. 

From Eq. 1)52(1 we can generalize that the overall phase depends on the flux of the magnetic field 
attached to the area of the loop. Given an arbitrary diagram we can initially calculate this flux to find 
out the overall phase. The typical nature of the overall phase function <i>(P, Q) is that P) is always 
unity. Similarly for all those graphs containing only two charged particle two-point functions, meeting 
each other at vertices P and Q, the overall phase function $(P, Q) = 1. But if the number of vertices 
in the loop, containing solely charged particles, is three or more then <f> is non trivial and it behaves like 
the Wilson loop as defined in Eq. 

5 Processes with charged particles as external legs in their Feyn- 
man diagrams 

Except the processes which we have discussed in the previous sections there are many other processes 
where there are charged particles in the external legs of the Feynman diagrams, as Compton scattering, 
electron electron-neutrino scattering, pair creation and annihilation and so on. In addition to the charged 
external particles all these processes may require virtual charged particles for their happening. In this 
cases the standard procedure is to solve the Dirac equation in presence of the background magnetic field. 
But in these cases the calculations loose manifest gauge invariance because we have to solve Eq. (jSJ 
choosing some form of the background gauge and that will fix the gauge in this case. To make the point 
clear we actually solve the Dirac equation in presence of a background magnetic field in this section. For 
stationary states, we can write the solution of the Dirac equation as: 



where 4> and x are 2-component objects. We use the Pauli-Dirac representation of the Dirac matrices. In 
this notation, we can write 




(52) 




(53) 



[E - m)<j) 
{E -I- m)x 



(T- (-iV-egAB)x, 
(T • (-iV - eqA^)(j), 



(54) 
(55) 
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where cr designates the Pauh matrices. EUminating we obtain 



cr ■ (— iV — eqA^) 



(56) 



Now we choose the background gauge field configuration as given in Eq. (|41|) with 6 = and with this 
choice Eq. (|56|l reduces to the form 



(57) 



Here is the diagonal Pauli matrix. With this choice of our background gauge field the positive energy 
solutions of the Dirac equation for an electron is given as: 

e-^P-^^C/,(2/,n,p^), (58) 

here denotes the space-time coordinate and p ■ = Ent — p ■ X^^ = Ent — PxX + PzZ. Here we have 
introduced the notation X for the spatial coordinates (in order to distinguish it from x, which is one of 
the components of X), and X^^ for the vector X with its y-component set equal to zero. In other words, 
p ■ X^ = PxX + PzZ, where Px and Pz denote the eigenvalues of momentum in the x and z directions.^ 
The energy En is given as: 



Ei 



771^ +pI + 2neB , 



(59) 



giving the relativistic form of Landau energy levels where n is a natural number. The form of Ug for 
s = +1, —1 are given by |12| . 









U+{y,n,p^) 



Pz 



E„ + m 
V2neB 



V 



In the above expressions 



and 



En 



( 



U-{y,n,p^) 



y/2neB 
En + m 

Pz 



E, 



(60) 



m 



/eB 



1/2 



where v are also natural numbers, and H^{S_) are Hermite polynomial functions of order v. 
The negative energy solutions are: 

e^P-''Ws{y,n,p,), 

where 



(61) 



(62) 



(63) 



En 



-in-m 



\ 



V2neB 



En + m 




( V2^ ^ \ 



V-iy,n,p,) 



En + m 

Pz 



En + m 




(64) 



^It is to be understood that whenever we write the spatial component of any vector with a lettered subscript, it would 
imply the corresponding contravariant component of the relevant 4-vector. 
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where ^ is obtained from ^ by changing the sign of the px-term. The calculation of the above solutions 
are given in the appendix at the end. 

The above solutions are obtained by solving the Dirac equation for a particular background gauge 
field configuration as specified in Eq. I|41|l with 6 = 0. If we designate the Dirac solutions in a background 
magnetic field by tpB{x) then the above solution is one of them but we can easily gauge transform 
this solution to obtain other solutions. All the gauge configurations we have specified in Eq. 1)41(1 . 
Eq. H42|) and Eq. (|43|l can be connected by smooth gauge transformations. Specifically if we take 6 = 
and c = then we can get the gauge fields in Eq. (|42|l from those of Eq. H41f) by the transformation 

— + + (9^A where A(x, y) — xyB. In this case "^p'si^) = i^B^x) g—^'^i^yB where iPb{x) is the solution 
for the gauge configuration as given in Eq. H41() and ip^ix) is the solution of the Dirac equation for the 
gauge configuration as given in Eq. H42I) for 6 = c = 0. In the above solutions the transverse momenta, 
Px and Py are just spurious degrees of freedom, they depend on the choice of the gauge. In the above 
solutions we will not find py as the gauge we chose to work with contained y, but had we started with the 
gauge fields as given in Eq. (|42|) then px should have been absent and if we had worked with the gauge 
fields as specified in Eq. (|43|l none of px or py should have appeared in our calculations. The arbitrariness 
of Px and Py reflects the fact that in presence of electromagnetic gauge fields Px and py are not the proper 
quantities to work with. But in what ever background gauge we choose to work with the expression of 
the energy as given in Eq. 159|l . which is a physically measurable gauge invariant quantity, remains the 
same. 

It can be a curious exercise to see if a pure gauge configuration like = Ag = Ag = and 
Ag = xk, where k is a constant, can produce any effect on the Dirac solutions. We saw in our previous 
discussion with the Schwinger two-point functions that a pure gauge configuration modifies its form 
but cannot affect loop calculations involving them. In the present case also it can be shown that pure 
gauge configurations although modifies the shape of the Dirac solutions but the free-particle solutions 
can be easily distinguished from the complete solutions. With the gauge which we have chosen we can 
proceed similarly to the case where we were solving the Dirac equation in a magnetic field, the equation 
corresponding to Eq. (|57|l in this case will be: 



2,^ - I ,^-eqxkX^ — - — 



(65) 



Looking at the above equation we can propose a possible solution as: 

= e'P-^^/(x), (66) 

where f{x) is a 2-component matrix which depends only on the ^-coordinate. In the present case 
p ■ = pyy + PzZ. Putting the above solution into Eq. (|^ we get the differential equation for f{x) 
which looks like: 

i—^eqxk\ f{x) = {E^-m'-pl~pl)f{x). (67) 

Now we can always assume a solution of the above equation of the form: 

f{x) = e"("+*'=«'=^)?7± , (68) 
where a is a constant and ri± are the standard two-component spinors of the form 

^^^^ 

Now if we put this solution into Eq. to get 













J. 




n 
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which predicts — m? —p^—pl — a? = and is the free-particle energy-momentum relation if we assume 
the constant a to be the momentum canonically conjugate to i.e. \i a = Px- So Eq. Hti6|l can now be 
written as: 




(71) 



which is similar to the free-particle solution up to the phase e*'^*'^^ . Using the expression of (j) we can find 
X from Eq. . It can be seen easily that the complete solution of the Dirac equation in the pure gauge 
which we chose is the free-particle solution times a phase. The above solution can also be interpreted in 
another way, all the solutions of the Dirac equation in presence of a pure gauge field are related to the 
free-particle solutions (solutions with = 0) by a gauge rotation. In the present case the rotation phase 
is simply e^'^^i'^^ . In presence of a magnetic field the Dirac solutions show that motion of the electrons 
transverse to the magnetic field directions gets quantized yielding discrit eigenvalues as n. This fact is 
crucial in making the B limit non trivial in the present circumstances. 

The solutions as given in Eq. and Eq. can be used to write the Dirac field in the usual 
sense and quantize the system in presence of a magnetic field ^1^]. Then using those fields one can 
define a two-point function of the fermions. But as the Dirac solutions which we obtained are solved for 
a particular gauge field configuration the calculations with them will not be manifestly gauge invariant 
but refiect the choice of the gauge. This situation is unlike the one which we faced when we were talking 
about the Schwinger two-point functions where gauge invariance of the loops can be manifestly proved. 
It may be tempting to write: 



where ^'^(a;) is the Dirac field in presence of a background magnetic field and is the appropriate 
vacuum state here, it is not the same as free-field vacuum as now the vacuum contains an infinite number 
of photons. iSBix,x') is the Schwinger two-point function as defined in Eq. (|24() . If we analyze the left 
and the right hand sides of the above equation we will see that such a relation cannot be correct. From 
the discussion presented in the end of section 13 it was observed that the Schwinger two-point function 
in the right is a function of A'^ and is defined for all gauge configurations including pure gauges and 
consequently it had dependence on the transverse momenta components p± and had a smooth i? — + 
limit. On the other hand 5's(x) which appears on the right hand side of Eq. H72|l is made up of Dirac 
solutions in a particular gauge which only follows Eq. ^ and as a result of this the p± components are 
spurious. There is no transformation which transforms Eq. H60|l or Eq. Ht)4|) to the pure gauge solutions 
as discussed in the previous paragraph. Precisely, the right hand side of Eq. H72|) is not a function of a 
general Aq{x) and this fact makes it different from the left hand side. As a consequence it is impossible 
to obtain the Schwinger two-point function by manipulating on the Dirac solutions which we obtain by 
solving Dirac equation in a magnetic field. 

From the discussions in this section we see that it is impossible to prove generally the gauge invariance 
of all those elementary particle processes which includes charged particles as external lines. The way of 
proving it must be particular and will depend upon the process under consideration. 

6 Conclusion 

In this article we tried to analyze the background gauge invariance of elementary particle processes in 
presence of a uniform background magnetic field. It was shown that the background gauge invariance 
of the calculations of certain class of processes requires a term in the fermion two-point functionss, in 
presence of a magnetic field, which must have proper covarince properties. In loops solely composed of 
charged particles the Schwinger phase will not contribute until the loop has three or more than three 
vertices. It was shown that pure gauge configurations will not alter any properties of the calculations and 
B ^ limit of the calculations yield the expected results. All the results stated are to one- loop but it is 
expected that the result is more general and can be prooved for higher loop cases also in a way similar 



iSb{x, x') - (i]|r{*s(x) ^B{x)}\n) , 



(72) 
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to the one presented in this article. Lastly we discuss all those particular processes which has one or 
more charged particles as external lines in their Feynman diagrams. In this cases it becomes impossible 
to furnish a general proof which will show that all these processes are background gauge invariant, as 
the calculations heavily rely on the effect of the particular gauge chosen to solve the Dirac equation in 
presence of a background magnetic field. The only way left to prove these processes to be background 
gauge invariant is to calculate a process in multiple gauges and show that the result is the same. 
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Appendix 

A Solution of the Dirac equation in a constant background mag- 
netic field 

In the standard Dirac-Pauli representation: 

where each block represents a 2 x 2 matrix. Noticing that the coordinates-ordinates x and z do not 
appear in Eq. H57I) except through the derivatives a possible solution of it can be: 

= e^^-^^/(2/), (74) 

where f{y) is a 2-component matrix which depends only on the y-coordinate, and possibly some momen- 
tum components, as we will see shortly. There will be two independent solutions for /(y), which can be 
taken, without any loss of generality, to be the eigenstates of with eigenvalues s = ±1. This means 
that we choose the two independent solutions in the form 

Since GzJs — sfs, the differential equations satisfied by Fg is 

^ - {eqBy + p,)^Fs + (E^ - ~ pi + eqBs)Fs = , (76) 
which is obtained from Eq. (|57|l . The solution is obtained by using the dimensionless variable 

(77) 



which transforms Eq. H76|l to the form 



F,^0, (78) 



where 



_ E^ -m^-pl + eqBs 
e\q\B ■ 

This is a special form of Hcrmitc's equation, and the solutions exist provided = 2v+l for v = 0,1,2, ■■ ■. 
This provides the energy eigenvalues 

= m'^ + pI + {2i^ + l)e\q\B - eqBs , (80) 
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and the solutions for are 

(81) 

where H^, are Hermite polynomials of order v, and N^, are normalizations which we take to be 

1/2 

(82) 

With our choice of normalization, the functions Ii, satisfy the completeness relation 




E UOiAC*) - ~ e*) = ^iy - y*) . (83) 

V 

where ^* is obtained by replacing y by in Eq. 1)77(1 . 

So far, q was arbitrary. We now specialize to the case of electrons, for which q = — 1. The solutions 
are then conveniently classified by the energy eigenvalues 

El=rr? ^'pl + lneB , (84) 

which is the rclativistic form of Landau energy levels. The solutions are two fold degenerate in general: 
for s — \^v — n— \ and for s = — 1, z/ = n. In the case of n = 0, only the second solution is available 
since v cannot be negative. The solutions can have positive or negative energies. We will denote the 
positive square root of the right side by i?„ . Representing the solution corresponding to this n-th Landau 
level by a superscript n, we can then write for the positive energy solutions, 

- ( '-f ) , ^ ( Zs) ) ■ («^' 

For n = 0, the solution /+ does not exist. We will consistently incorporate this fact by defining 

/-i(y)-o, (86) 

in addition to the definition of /„ in Eq. ((81|l for non-negative integers n. 

The solutions in Eq. H85|l determine the upper components of the spinors through Eq. (|74|l . The lower 
components, denoted by x earlier, can be solved using Eq. ifS^ . 
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